Mon. Not. R. Astron. Soc. 000, [T]-?? (2005) 



Printed 5 February 2008 



(MN IATeX style file v2.2) 



Rotational velocities of the giants in symbiotic stars: 
II. Are S-type symbiotics synchronized? * 

R. K. Zamanov^'^ f, M. F. Bode^, C. H. F. Melo3'^ R. Bachev\ 
A. Gomboc^'S, I. K. Stateva\ J. M. Porter^, J. Pritchard^ 

^ Institute of Astronomy, Bulgarian Academy of Sciences, 72 Tsarighradsko Shousse Blvd., 1784 Sofia, Bulgaria 

^ Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Birkenhead, CH4I ILD, UK 

^ European Southern Observatory, Casilla 19001, Santiago 19, Chile 

* Departamento de Astronomia, Universidad de Chile, Casilla 36-D, Santiago, Chile 

^ Department of Physics, University of Ljubljana, Jadranska 19, 6100 Ljubljana, Slovenia 



Accepted . Received 2007 May 17; in original form 2006 March 17 



ABSTRACT 

We have measured the projected rotational velocities (u sini) of the mass donors 
for 29 S-type symbiotic stars using high resolution spectroscopic observations and the 
cross-correlation function (CCF) method. The results of the CCF have been controlled 
with synthetic spectra. The typical rotational velocity of the K and M giants in S- 
type symbiotics appeared to be 4.5 <v sini< 11.7 kms~^. In a sub-sample of 16 
S-type symbiotic stars (with known orbital periods and well measured v sini) 15 have 
deviations from synchronization less then the 3-(t level. This means that wc did not 
find evidence for a statistically significant deviation from the synchronization for any 
of these 15 objects. The deviation from synchronization is statistically significant (at 
confidence level >99%) only for the recurrent nova RS Oph. 

For 22 S-type symbiotics we give clues as to what their orbital periods could be. 
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1 INTRODUCTION 

The Symbiotic stars (SSs - thought to comprise a white 
dwarf (WD) accreting from a cool giant or Mira) represent 
the extremum of the interacting binary star classification. 
They offer a laboratory in which to study such important 
processes as (i) mass loss from cool giants and the formation 
of Planetary Nebulae; (ii) accretion onto compact objects, 
(iii) photoionisation and radiative transfer in gaseous neb- 
ulae, and (iv) nonrelativistic jets and bipolar outflows (e.g. 
Kenyon 1986; Corradi et al. 2003). 

On the basis of their IR properties, SSs have been clas- 
sified into stellar continuum (S) and dusty (D or D') types 
(Allen 1982). The D-type systems contain Mira variables as 
mass donors. The D'-type are characterized by an earlier 
spectral type (F-K) of the cool component and lower dust 
temperatures. All mass donors in D'-type systems appeared 
to be very fast rotators (see Paper I, Zamanov et al. 2006). 

Our aims here are: (1) to measure the projected rota- 
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f e-mail: rkz@astro.bas.bg; mfb@astro.livjm.ac.uk; 

andrej a.gomboc@fmf . uni-lj . si 



tional velocities {v sini) and the rotational periods {Prot) of 
the giants in a number of southern S-type SSs, using a cross 
correlation function (CCF) approach; (2) to check whether 
their rotation is synchronized with the orbital period; (3) 
to provide pointers to the determination of binary periods 
(assuming co-rotation). 

This is the second in a series of papers exploring the 
rotation velocities of the mass donating (cool) components 
of SSs. 



2 OBSERVATIONS 

We have observed 30 objects from the Belczynski et al.(2000) 
SS catalogue and have observed all S-type SSs from the cata- 
logue with 12'' < RA < 24'', declination < 2*^, and catalogue 
magnitude brighter than V < 12.5. 

The observations have been performed with FEROS 
at the 2.2m telescope (ESO, La SUla). FEROS is a fibre- 
fed echelle spectrograph, providing a high resolution of 
A/ A A =48000, a wide wavelength coverage from about 
4000 A to 8900 A in one exposure and a high throughput 
(Kaufer et al. 1999). The 39 orders of the echelle spectrum 
are registered with a 2k x 4k EEV CCD 
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Table 1. Journal of observations. In the table are given as follows: 
the name of the object, date of observation (YYYY-MM-DD), 
the modified Julian Date (JD - 2400000.5) of the start of the 
observation, signal to noise ratio (S/N) around A 8500 A. 



object date-obs MJD-OBS Exp S/N 

[sec] 

40 
45 
45 
40 
70 
45 
80 
40 
70 
50 
45 
40 
75 
50 
35 
40 
55 
70 
45 
50 
80 
65 
50 
45 
40 
90 
40 
50 
50 
50 



The present data have been collected from April to 
September 2004. Table [T] gives a log of the observations. 
All spectra are reduced using the dedicated FEROS data 
reduction software implemented in the ESO-MIDAS system 
(www.ls.eso.org/lasilla/sciops/2p2/E2p2M/FEROS/DRS/). 
A few examples of our spectra are given in the Appendix 
(Figll. 



3 « sini MEASUREMENT TECHNIQUES 
3.1 CCF method 

The projected rotational velocities have been derived by 
cross-correlating the observed spectra with KG numerical 
masks yielding a cross-correlation function whose width 
(cTobs) is related to broadening mechanisms such as stellar 
rotation and turbulence. 

The emission lines do have an effect in the CCF and 
they must be cleaned. They were cut off by fitting a con- 
tinuum and replacing the emission lines by the value of the 
fit. Note that the exact location of the continuum is not 
important. 

The numerical KO-mask was constructed from a KOIII 



synthetic spectrum in the region between AA 5000-7000 A 
following the procedure described in Baranne et al. (1979) 
In the CORAVEL-type cross-correlation a binary mask is 
used as template instead a real spectrum. This binary 
(or CORAVEL-type) mask has been used in many differ- 
ent data reduction software for cross-correlation (ELODIE, 
CORALIE, and now HARPS). The KO-mask CCFs for the 
SSs observed here are plotted in the Appendix (FiglSj and 
Gobs are given in Tabled 

In order to use the observed width of the CCF {(Jobs) 
as an estimate of v sini one needs to subtract the amount 
of broadening contributing to aobs unrelated to the stellar 
rotation (convection, instrumental profile, etc), i.e. ctq. For 
FEROS spectra (see Melo et al. 2001, 2003): 

u sini = 1.9 \/ (^l^^ - (y\ kms"^ (1) 

More details of the cross-correlation procedure are given 
in Melo et al. (2001), and also cro is calibrated as a function 
of the (B — y) for FEROS spectra and for stars with 0.6 < 
{B -V) < 1.2. For giants with (B -V) < 1.2 in Table [2] 
as well as in Paper I, the Melo et al.(2001) calibration has 
been adopted. However, the stars in our sample have (B — V) 
around 1.5 which is beyond the range of the calibration of 
Melo et al.(2001). For (B-V)> 1.2, we wiU adopt constant 
ao = 4.5 kms~^. 

This value has been adopted as a result of: (i) CCF 
measurements of few objects with known v sin i; (ii) bearing 
in mind the results of Delfosse et al. (1998). Using a similar 
template and method but lower resolution (which increases 
(To), Delfosse et al. (1998) have shown that ao does vary over 
a range of 0.8 < {R~ I) < 1.5 (which corresponds to a spec- 
tral type from ~M0 to ~M6), decreasing from ~ 5.1 kms~^ 
to 4.7 kms~^; (iii) bearing in mind the calibration of CCF 
and FEROS spectra undertaken in Melo et al. (2001). 

The errors on CCF v sin i measurements are dominated 
by systematic effects rather than by photon noise. The error 
on V sini comes from two main sources: uncertainties on the 
values of aobs and ao. According to error estimate carried out 
by Melo et al.(2001), stars with (B-V) < 1.2 have an error 
on V sin i less than 1.2 kms~^. For stars with {B — V) > 1.2 
we adopt a constant cro, which leads to an increase in our 
errors. For such stars a conservative error of 1.5 kms^^ is 
assigned. For objects with v sin 15 kms~^ the error on 
V sini is ±10%. 

In the case of the rapidly rotating Hen 3-1674 the CCF 
rotation was extracted by a slightly different procedure as 
described for the fast rotators in Paper I. For V3804 Sgr we 
did not get a meaningful CCF due to the numerous emis- 
sion lines in the spectrum at the time of our observations. 
From measurements of the Fel 8689 line we obtain a rough 
estimate of FWHM(FeI 8689)^ 0.58 ± 0.15 A similar to the 
width of this line in AR Pav and EG Ser and corresponding 
to V sin i« 9 ± 3 km s~^ . 



3.2 FWHM method 

Besides the cross-correlation procedure another method us- 
ing full width at half maximum (FWHM) of spectral lines of 
observed and synthetic spectra was also applied. This proce- 
dure is similar to that described in Fekel (1997). The spec- 
tra of K5III (Te//=3950 K, logg = 1.5), MOIII (Te/7=3985 
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Hen 2-374 


2004- 


-06- 


-08 


53164, 


.3099 


1200 


Hen 3-1213 


2004- 


-04- 


-12 


53107, 


.3049 


1200 


Hen 3-1341 


2004- 


-04- 


-12 


53107, 


.3354 


1200 


Hen 3-1674 


2004- 


-06- 


-07 


53163, 


.3103 


1200 


Hen 3-1761 


2004- 


-06- 


-04 


53160, 


.2095 


600 


Hen 3-863 


2004- 


-04- 


-11 


53106, 


.3067 


1200 


MWC 960 


2004- 


-06- 


-08 


53164, 


.2738 


1200 


PN Ap 1-9 


2004- 


-06- 


-07 


53163, 


.1986 


1200 


RS Oph 


2004- 


-04- 


-11 


53106, 


.3849 


600 


RW Hya 


2004- 


-04- 


-12 


53107, 


.2285 


600 


SS73 129 


2004- 


-06- 


-29 


53185, 


.0660 


1200 


SS73 141 


2004- 


-06- 


-29 


53185, 


.1058 


1200 


V2506 Sgr 


2004- 


-06- 


-07 


53163, 


.2393 


1200 


V2756 Sgr 


2004- 


-06- 


-08 


53164, 


.3968 


600 


V2905 Sgr 


2004- 


-06- 


-28 


53184, 


.0858 


1200 


V3804 Sgr 


2004- 


-06- 


-27 


53183, 


.9951 


1200 


V4018 Sgr 


2004- 


-08- 


-31 


53248, 


.1800 


1200 


V4074 Sgr 


2004- 


-06- 


-07 


53163, 


.2737 


600 


V919 Sgr 


2004- 


-06- 


-03 


53159, 


.3108 


1200 
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Table 2. Rotational velocities of the red giants in symbiotic stars (measured in this paper). The spectral types are from the catalogue of 
Belczyriski et al.(2000), (B-V)o is the intrinsic colour (from Schmidt-Kaler 1982) for the corresponding spectral type adopting luminosity 
class III in all cases, Kg is the adopted radius of the giant [if not otherwise indicated, this is the average radius for the corresponding 
spectral type taken from Table 7 of van Belle et al.(1999)]. Port is the orbital period (see Sect[6]l. v sini(FWHM) is the projected 
rotational velocity of the cool giant as measured with FWHM method (see Sect. 3. 2). <Toi,a is the observed width of the CCF. v sini is 
our measurement based on CCF method (see Sect. 3.1), if other measurements of v sini exist, they are also given in the next column. 
For objects with unknown orbit we give an estimation of the upper limit of the orbital period (i.e. we expect Port ^Puii see Sect|9ll. 
The upper part of this table contains objects observed with FEROS, the lower part - v sini values are taken from the literature. 



object 


Cool 


(B-V)o 


Rg 


Porb 


V sini 


<^obs 


V sini 


other 














FWHM 


CCF 


CCF 








Sp. 




R© 


[days] 


km s~^ 


[kms"-'-] 


[km s^-*-] 


[km s^-*-] 


[days] 


AR Pav 


M5 


1.59 


139.6 


604.5 


8.8±2 


6.208 


8.1±1.5 


11±2'= 


— 


AS 255 


K3 


1.26 


20.5 




9.7±1 


6.417 


8.7±1.5 




119 


AS 276 


M4.5 


1.60 


123.0 




7.7±2 


6.382 


8.6±1.5 




724 


AS 289 


M3.5 


1.62 


89.0 


451 


9.5±1 


8.400 


13.5±1.5 


5.7±1* 


— 


AS 316 


M4 


1.62 


105.5 




9.6±1 


6.843 


9.8±1.5 




545 


AS 327 


M3 


1.62 


71.5 




7.7±1 


5.847 


7.1±1.5 




510 


BD-21°3873 


K2 


1.16 


20.8 


281.6 


6.6±1 


5.352 


4.6±1.2 


5.4±0.7^ 


— 


CD-36°8436 


M5.5 


1.58 


144.0 




6.6±1.5 


6.410 


8.7±1.5 




840 


CD-43° 14304 


K5 


1.51 


38.8 


1448 


7.2±1 


5.869 


7.2±1.5 


<39 


— 


FG Ser 


M5 


1.59 


139.6 


650 


9.8±1 


6.756 


9.6±1.5 


8±l'', 7±l'= 


— 


HD 319167 


M3 


1.62 


71.5 




7.4±2 


6.041 


7.7±1.5 




472 


Hen 2-374 


M5.5 


1.59 


144.0 




8.4±1.5 


5.712 


6.7±1.5 




1091 


Hen 3-1213 


M2/K4 


1.54 


50:? 




9.1±1.5 


7.236 


10.8±1.5 




235 


Hen 3-1341 


M2 


1.61 


57.8 




7.5±1.5 


6.423 


8.7±1.5 




336 


Hen 3-1674 


M5 


1.59 


139.6 




56.0±5 




52.0±5.2 




135 


Hen 3-1761 


M4 


1.62 


105.5 




9.3±2 


6.475 


8.8±1.5 




603 


Hen 3-863 


K4 


1.43 


45.0 




7.9±1 


5.811 


7.0 ±1.5 




326 


MWC 960 


K9 


1.55 


39.0" 




9.0±1 


6.436 


8.7±1.5 




225 


PN Ap 1-9 


K4 


1.43 


45.0 




8.2±1 


5.792 


6.9±1.5 




329 


RS Oph 


MO 


1.56 


59.1 


455.7 


13.8±1.5 


7.624 


11.7±1.5 






RW Hya 


M2 


1.61 


57.8 


370.2 


6.2±1 


5.837 


7.1±1.5 


5.0±1'= 




SS73 129 


MO 


1.56 


59.1 




8.9±1 


6.162 


8. Oil. 5 




374 


SS73 141 


M5 


1.59 


139.6 




7.9±1 


6.064 


7.7±1.5 




914 


V2506 Sgr 


M5.5 


1.59 


144.0 




7.8±1 


6.284 


8.3±1.5 




874 


V2756 Sgr 


M3 


1.62 


71.5 


243:? 


4.2±1.5 


4.942 


3.9±1.5 




932 


V2905 Sgr 


M5 


1.59 


139.6 




6.6±1 


5.708 


6.7±1.5 




1059 


V3804 Sgr 


M5 


1.59 


139.6 




9.0±3 










V4018 Sgr 


M4 


1.62 


105.5 




7.2±1.5 


5.345 


5.5±1.5 




974 


V4074 Sgr 


M4 


1.62 


105.5 




4.2±1.5 


4.871 


3.5±1.5 




1508 


V919 Sgr 


M2 


1.61 


57.8 




7.4±1 


6.560 


9.1±1.5 




322 



SY Mus 


M5 


139.6 


624.5 


7±1P 


AG Dra 


K2Ib or II 


30-40'' 


554 


3.6±l'=, 5.9±1'' 


BX Mon 


M5III 


139.6 


1401 


6.8±l'= 


TX CVn 


K5III 


38.8 


199 


8.9±l''' 


T CrB 


M4IIIelips. 


66±11= 


227.57 


5.4±l'= 


V443 Her 


M5.5III 


144 


594 


4.5±l'= 


CI Cyg 


M5II or M5.5 


189-236"* 


854.5 


10.4±1* 


AG Peg 


M3III 


71.5 


817.4 


4.5±l'= 


V1329 Cyg 


M6 


147.9 


956.5 


7.0±2*= 


BF Cyg 


M5III 


75-280'' 


756.8 


4.5±2'= 



"-from the fit in van Belle et al.(1999); ''Tomov et al.(2000); "^Belczynski & Mikolajewska(1998) ''Kenyon et al.(1995); '=Schild et 
al.(2001); ■''Smith et al.(1997); ^Schmid ct al.(1998); '^Miirset et al.(2000); '"Fekel et al.(2003); ^Schmutz et al.(1994); ^de Medeiros & 

Mayor (1999). ""see Sectlg^] 



K, logg = 1.2), and M5III (Te//=3424 K, logg = 0.5), 
stars have been synthesized by using the code SYNSPEC 
(Hubeny, Lanz & Jeffery 1994) in the spectral region A 8750- 
8850 A. 

LTE model atmospheres were extracted from Kurucz's 



grid (1993). The VALD atomic line database (Kupka et al. 
1999) was used to create a line list for spectrum synthesis. 
The value of 3 kms^^ was adopted for the microturbulent 
velocity. A grid of synthetic spectra for projected rotational 
velocities from kms~^ to 60 kms~^ was calculated. The 
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5.0 10.0 15.0 

V sin i (CCF) [km s-i] 



Figure 1. Comparison between the CCF and FWHM 
measurements of v sini. The straight line indicates 
V sinj(CCF)=i; sinj(FWHM). 



FWHM of a dozen observed spectral lines has been mea- 
sured and compared to the FWHM of spectral lines from 
the synthetic spectra. The results are given in the sixth col- 
umn of Table El 

The comparison between the CCF and FHWM methods 
(see Figd]) shows that the v sini measurements agree well, 
with typical difference ±1 — 2 kms~^, as expected from the 
measurement errors. The only exception is AS 289 (see also 
Sect. 6.1). 

Because the CCF method uses many more lines and a 
better defined mathematical procedure, in the analysis we 
will use V sin i derived with the CCF method. 



4 PARAMETERS OF SYMBIOTIC STARS 
4.1 Synchronization in SSs 

The physics of tidal synchronization for stars with convec- 
tive envelopes has been analyzed several times (e.g. Zahn 
1977, and see the discussion in Chapter 8 of Tassoul 2000) . 
There are some differences in the analysis of different au- 
thors, leading to varying synchronization timescales. Here, 
we use the estimate from Zahn (1977, 1989). The synchro- 
nization timescale in terms of the period is 



800 



v3m,^(|| + i) 



orb 



years 



(2) 



where Mg and M2 are the masses of the giant and white 
dwarf respectively in Solar units, and Rg and Lg are the 
radius and luminosity of the giant, also in Solar units. The 
orbital period Porb is measured in days. 

In Table [3] we calculate this timescale for a few rep- 
resentative cases. We point out that these values serve as 



Table 3. The synchronization time, Tsyn, calculated following 
Eq.l for a few representative cases. The typical values for Lg and 
Rg of K5 III and M4III stars are adopted following van Belle et 
al.(1999). 



Mass donor 




Mg 


Rg 




M2 




Sp. type 


[Lq] 


[Mq] 


[Rq] 




[Mq] 


[yr] 


K5 III 


350 


2 


40 


200 d 


1.0 


7000 










1500 d 


1.0 


2.10^^ 


M4 III 


1380 


1 


105 


500 d 


0.5 


140 






2 




500 d 


1.0 


700 






3 




500 d 


1.0 


3200 






2 




1500 d 


1.0 


60000 



an estimation. From the Rg dependence of Tsyn in Eq. 1, 
one can expect large uncertainties in the synchronization 
timescale as a result of the uncertainties in the red giant ra- 
dius. There is also a P^rb dependence, which can introduce 
large uncertainties, but Porb is much more accurately known 
than Rg. 

As can be seen in Table |3l for S-type SSs depending 
on the binary parameters, Tsyn can be as short as < 10'' yr 
and longer than lO'^ years. Short Tsyn can be expected when 
the radius of the mass donor is a significant fraction of the 
orbital separation. 

The lifetime of the symbiotic phase for a red giant or 
AGB star is around lO^year (see e.g. Yungleson et al. 1995), 
which is comparable with Tsyn- 



4.2 Inclination 

In our calculations, we assume that the rotational axis of 
the red giant is perpendicular to the orbital plane (see Fekel 
1981; Hale 1994; Stawikowski 1994). To calculate the rota- 
tional periods of the mass donors (Prot) we need to know 
the inclination of the orbit (i) and Rg. 

For eclipsing binaries we can adopt inclination i = 70'^ — 
110", which produces small errors in sini. For such cases we 
will assume that sini is in the range 0.94-1.00. For other 
cases we use results from spectropolarimetric observations 
or radial velocity measurements. 



4.3 Radius of the cool component 

For a few objects, the radii of the cool components are de- 
rived from model calculations. If the radius of the giant (Rg) 
is not known, we will use the average radius for the corre- 
sponding spectral type taken from van Belle et al.(1999), 
always adopting luminosity class HI and error ±5%. 

For the M giants these values are similar to the stellar 
radii of M giants in the Hipparcos catalogue as calculated 
by Dumm & Schild (1998). Dumm & Schild (1998) have 
also shown that the radius of the M giants depends on the 
mass. For 9 objects included in our sample, the masses of 
the M giants are listed in Mikolajewska(2003). In Table HI 
we compare the values adopted here (following van Belle et 
al. 1999) with the results of Dumm & Schild (1998). 

As visible in most cases, the adopted radii (following 
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Table 4. Stellar radii of M giants. In column 3 are given Kg 
(following van Belle et al. 1999), in column 4 - mass of the M 
giants (from Mikolajewska 2003), in column 5 - the range of Kg 
for the corresponding mass and spectral type (following Dumm 
& Schild 1998). 



Object 


Mass donor 


Kg 


Mg 

[M©] 


Kg 






[R©] 


[R©] 


1 


2 


3 


4 


5 


FG Ser 


M5III 


140 


1.7 


110-150 


AR Pav 


M5III 


140 


2.5 


90-240 


AG Peg 


M3III 


72 


>1.8 


60-90 


BX Mon 


M5III 


140 


3-3.7 


150-200 


SY Mus 


M5III 


140 


1.3 


60-100 


RW Hya 


M2III 


58 


1.6 


40-80 


AR Pav 


M5III 


140 


2.5-3 


130-160 


BF Cyg 


M5III 


140 


1.8 


110-140 


AG Peg 


M3III 


72 


^1-8 


60-90 



van Belle et al. 1999) are in agreement with the mass- 
dependent radii of Dumm & Schild (1998). There are 2 cases 
where the adopted values are outside of this range: BX Mon 
and SY Mus (see SectO. 



5 RESULTS 

Our measurements of v sini together with data collected 
from the literature, are summarized in Table [21 Our sample 
thus contains 39 objects (29 with ii sini measured by us and 
10 taken from the literature). This sample should have no 
biases in the rotational speed of the cool giant, even though 
the sample is flux limited. 

In our sample of S-type SSs, the projected rotational 
velocities of the mass donors are from 3.5 kms~^ up to 
52 kms~^. Including our data and the data from the litera- 
ture, we obtain for v sini: mean=8.7 kms~^, median= 7.7 
kms~^, standard deviation^ 7.4 kms~^. 

If we exclude the 2 slowest and the 2 fastest rotators we 
get: mean v sini= 7.6±1.8 kms~^. In effect, 90% of the mass 
donors in S-type SSs with measured rotation have v sini in 
the interval 4.5 sini^ 11.7 kms~^. For the southern S- 
type symbiotics (29 objects, flux limited sample, FEROS ob- 
servations) the values are similar: median v sini=8.1 kms~^, 
90% of the objects in the interval 4 sini^ 11.7 kms~^. 



6 S-TYPE SYMBIOTICS WITH KNOWN 
ORBITAL PERIODS 

In our sample, there are 17 objects with known orbital peri- 
ods and 1 where this is inferred (V2756 Sgr) . In this section 
we compare the orbital periods with the rotational periods 
of the mass donors object by object. The upper-lower limits 
of Prot are calculated from: 

_ 27r(_Rg ±ei)(sini±e2) 

J rot — : — : , [o) 

V sm I =F 63 

where ei, 62, and 63 are the corresponding errors in Kg, sini, 
and V sini respectively. 



In the statistical analysis (Sect|8}, for Prot we will use 
the average value between these upper-lower limits and will 
consider them as corresponding to ±la error. 

6.1 Objects with known orbital periods observed 
by us 

AR Pav is an eclipsing binary with Por6 ~ 604.5 days 
(Bruch et al. 1994). The giant's radius derived from the 
eclipse is Kg = 137 ± 20 Rq (Quiroga et al. 2002), which is 
very close to the value adopted in Table [2] from the spectral 
type of the cool giant and to Rg = 139 ± 10 -Rq (Skopal 
2005). We calculate 649 <Prot < 1120 days. 

AS 289 (V343 Ser) has orbital period Por6 = 450.5 ± 

2.2 d, e = 0.135 ± 0.046, and orbital inclination of 17° - 23° 
(Fekel et al. 2001). Kenyon & Fernandez-Castro (1987) gave 
spectral type M3.9III, but more recently M3.5 is assigned 
from Miirset & Schmid (1999). Assuming it is a normal 
M3.5III star, and using our CCF value for v sini, we ob- 
tain 159 < Prot < 341 days. 

In a binary with an eccentric orbit the synchronization 
is reached at a value of Prot <Por6 (pseudosynchronization, 
see Hut 1981). Following Eq. 42 and Eq. 43 of Hut (1981), 
in AS 289 the pseudosynchronization is expected at Prot ~ 

0.90 Porb. 

It is worth noting that Fekel et al.(2003) give a value of 
V sini= 5.7 ±1 kms~^ and our FWHM method give v sini= 
9.5 ± 1 kms^^. An independent measurement of v sini with 
better spectroscopic resolution will be valuable. 

BD-21°3873 has P^rb = 281.6 ± 1.2 days, sin i = 0.87 
(Smith et al. 1997). We derive 147 < Prot < 293 days. The 
object is thus synchronized within the measurement errors. 

RW Hya is an eclipsing binary with Porb = 370.4±0.8 
days (Schild et al. 1996). Supposing i = 70" - 90°, we cal- 
culate that Prot is in the interval 301 - 552 days. The object 
is thus synchronized within the measurement errors. 

FG Ser (AS 296) is an eclipsing binary with Porb = 
650 ± 5 days (Miirset et al. 2000). Supposing i = 70° - 
90°, we calculate Prot ~ 563 — 918 days. The object is thus 
synchronized within the measurement errors. 



6.2 Objects with known periods, not observed by 
us 

Fekel, Hinkle & Joyce (2004) published values of v sini for 
13 S-type SSs. In this section, if not otherwise stated, we 
use their measurements. 

SY Mus is an eclipsing binary with i = 95°. 8 ± 1.7° 
(Harries & Howarth 1996, 2000), Porb = 624.5 days and 
V sini = 7±1 kms~^ (Pereira et al.l995; Schmutz et al. 
1994; Kenyon & Mikolajewska 1995). We calculate Prot ~ 
830-1230 days. [If we use Kg x 80 Rq (see Sect|43I, we get 
Prot ~ 480 — 705 days, a range which includes the orbital 
period.] 

AG Dra is a yellow SS with spectroscopic Porb ~ 
548.5±2 d (Friedjung et al. 2003; Fekel et al. 2000b), and 
i r; 30° - 45° (Mikolajewska et al.l995). The cool com- 
ponent is a K2 lb or II (Zhu et al. 1999). There are two 
measurements of v sini: 5.9±1.0 kms"^ (de Medeiros & 
Mayor 1999) and 3.6±1.0 kms~^ (Fekel et al. 2004). We 
adopt V sini= 4.8 ± 1.0 kms'V Skopal (2005) calculated 
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Figure 2. The rotational period of the red giant (Prot) versus the orbital period {Porb) of the 17 objects in our sample. The solid line 
corresponds to Prot =Porb- Most objects are close to this line, which indicates that they are synchronized. There are 2 objects which 
deviate considerably from that rule (RS Oph and CD-43° 14304). 



Rg = 33 ± 11 Rq. Following Tomov et al.(2000) and the 
references therein we can adopt Rg = 30 — 40 Rq . We then 
estimate Prot ~ 130 — 485 days. 

TX CVn has Vorb = 199±3 d, e = 0.16±0.06, with the 
cool component a normal K5III star, and inclination 20° < 
i < 70° (Kenyon & Garcia 1989). The luminosity class could 
also be II or lb (Zhu et al. 1999). We calculate P^ot « 64 - 
242 days. Pseudosynchronization (see Hut 1981) is expected 
at Prot ~ 0.87Port). The object is thus (pseudo) synchronized 
within the measurement errors. 

AG Peg has Porb = 818.2 ± 1.6 days and e = 0.110 ± 
0.039 (Fekel et al. 2000a), normal M3III and i ~ 40° - 60° 
(Kenyon et al. 1993). We calculate Prot ~ 400 - 933 days. 
The object is thus synchronized within the measurement 
errors. 

V1329 Cyg has Pori, = 963.1 ± 9.8 days (Chochol & 
Wilson 2001) and inclination i = 86" ±2" (Schild & Schmid 



1997). We calculate Prot ~ 734 - 1571 days. The object is 
thus synchronized within the measurement errors. 

T CrB has Porb = 227.57 d (Fekel et al. 2000a), 
and Rg = 66 ± 11 i?Q (Belczynski & Mikolajewska 1998). 
[It should be noted that Skopal (2005) gives Rg = 75 ± 
12(d/960pc), and that a normal M4III would have Rg ~ 
105. 5R0.] The inclination of the system is i ~ 65" — 70" 
(Stanishev et al. 2004). 

We calculate Prot ~ 391 — 832 days, which differs from 
the expectations for synchronization. This result is unex- 
pected, however it appears that this deviation is not statis- 
tically significant (see Table[S}. From Eq.l we get Tsyn < 100 
yr. The red giant in T CrB is ellipsoidally shaped (Yudin & 
Munari 1993). This means that the object has to be syn- 
chronized. 

BX Mon The orbital parameters of this eclipsing 
system arc: 

Porb = 1401 days, eccentricity e — 0.49 (Dumm et al. 1998); 
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Port =1259 d, e = 0.44 (Fekel et al. 2000a). 

We calculate Prot ~ 800 — 1278 days [if we use Rg « 
170 R© (see SectHlil, we get Prot ~ 970 - 1550 days]. 
BX Mon is the only SS with considerable orbital eccentricity 
in our sample. SSs with Port > 800 days tends to eccentric 
orbits (Fekel et al. 2007). In a binary with e = 0.44, the pseu- 
dosynchronization is expected at about Prot ~ 0.46 Porb (see 
Hut 1981). 

The time scale for synchronization in SSs is ~10 times 
shorter then the circularization time (see Schmutz et al. 
1994). All this implies that in this system the red giant 
is more or less (pseudo)synchronized, but the orbit is not 
circularized yet. It has therefore to be in the process of cir- 
cularization. 

V443 Her is not eclipsing, with Port = 599.4±2.1 days 
(Fekel et al. 2000b), viewed at an inclination i ~ 30° (Do- 
brzycka et al. 1993). We calculate Prot = 528 - 1245 days. 
The red giant is thus synchronized within the measurement 
errors. 

CI Cyg is an eclipsing binary with orbital period 855.6 
days and orbital separation a — 2.2 au (Kenyon et al. 1995). 
Kenyon et al.(1995) have shown that the mass donor is 
an M5II asymptotic branch giant, filling its tidal surface. 
They calculated from the eclipse that Kg/a ~ 0.4 — 0.5 
which means that Rg « 189 — 236 Kq. We calculate Prot ~ 
780 — 1270 days. The object is thus synchronized within the 
measurement errors. 

BF Cyg is an eclipsing binary with Port = 757.3 days 
and inclination i ~ 70 — 90^^ (Pucinskas 1970, Skopal et 
al. 1997). If the cool component is a normal M5III giant 
we expect Rg = 139.6 ± 5% R©. Mikolajewska et al.(1989) 
suggested Rg ~ 75 R©, however Skopal et al. (1997) give 
Rg = 260±20 Rq. Using v sini= 4.5±2 kms"'^, we calculate 
Prot ~ 560 - 5660 days. 

There are signs that the red giant is ellipsoidally shaped 
(Yudin et al. 2005). Tsyn for such an object will be short, 
from Eq.l we get Tsyn < 7000 yr. It means that the red 
giant has to be tidally locked and most probably the Rg is 
closer to a value ~ 75 Kq. 



7 QUESTIONABLE OBJECTS 

RS Oph has orbital period Porb = 455.72±0.83 days (Fekel 
et al. 2000a). The inclination is about 30-40° (Dobrzycka 
& Kenyonl994; Dobrzycka et al. 1996). This will give a 
108 <Prof < 197 days, i.e. 2-3 times less than the orbital 
period. 

RS Oph is a peculiar SS exhibiting different types of 
activity - recurrent nova eruptions, jet or blob ejections, 
flickering (see Bode et al. 2006 and references therein). 

This is one of the two objects in our sample whose devi- 
ation from synchronization is statistically significant (at con- 
fidence level >99%, see Table [5]). The red giant in RS Oph 
seems to rotate faster than the orbital period, which means 
that it has to be in a process of deceleration (the expected 
Tsyn «; 5.10* yr). 

While there are no doubts about Porb, any of the other 
parameters {v sini, inclination, red giant radius) have to be 
checked with independent measurements. 

It is noteworthy that our experiments to measure v sin i 



with CCF and spectra from different epochs showed that 

V sini could even be higher (up to 14.5±1.5 kms"'^). 

CD-43°14304 Schmid et al.(1998) reported a circular 
orbit, Porb = 1448±100 days. They argued that the presence 
of phase-dependent Ha variations is attributable to occulta- 
tion effects, suggesting a relatively high inclination, i > 45''. 
Spectropolarimetry (Harries & Howarth 2000) gives two pos- 
sible values for the inclination of the system i = 57° ± 5° 
or i = 122'' ± 48". The second corresponds to an improba- 
bly large eccentricity for the orbit. Therefore we can accept 
sini = 0.79 — 0.88. With our value of v sini we calculate 
Prot ~ 170 - 321 days. 

However, Schmid et al.(1998) found no evidence for ro- 
tational broadening in the cool-giant spectrum and were able 
to place an upper limit on v sini< 3 kms~^. We measured 

V sini= 7.2±1.5 kms~^ the CCF looks good (see Fig.O and 
FWHM method gives a similar result. If our CCF result is 
wrong and that of Schmid et al.(1998) is a correct one, then 
the object is probably synchronized. To avoid confusion we 
will exclude this star from the analysis, but an independent 
check of 11 sini would be valuable. 

V2756 Sgr (Hen 2-370, SS73 145). The parameters of 
the system are not well known. We calculate Prot < 1595 
days. A photometric period of 243 days is supposed in Hof- 
fleit (1970). This photometric period is not confirmed with 
radial velocity measurements to be that of the orbit. If this 
is the orbital period and the red giant is tidally locked, then 
the inclination of the system would be about i ~ 15". 

Hen 3-1674 is the fastest rotator in our sample. The 
catalogues indicate that the mass donor is probably M5III 
star (see Table[2]and the references therein). A normal M5III 
star would have Rg « 90—160 Kq (van Belle et al 1999) and 
Mg 1 — 3 Mq, which means a break-up velocity of 30-60 
km s~^ . If our measurements and the adopted parameters of 
the system are correct then Hen 3-1674 rotates close to its 
critical velocity and has rotation similar to that of D'-type 
SSs (see Paper I). 



8 ARE THE MASS DONORS IN S-TYPE 
SYMBIOTICS SYNCHRONIZED 
(CO-ROTATING) ? 

Fig. [2] shows the rotational period versus the orbital period 
of the 17 objects in our sample, with a straight line indicat- 
ing the co-rotation (i.e. Prof=Porb). Most objects are close to 
this line, which suggests that they are synchronized. In Ta- 
ble [S] are given the individual deviations as well as the corre- 
sponding probability that the deviation is random. 9 objects 
are synchronized within the measurement errors (1-cr level). 
4 objects have deviations between 1 and 2-a. Generally, 15 
out of 17 are within the 3-a level. The two objects that are 
outside of the 3-a level are RS Oph and CD-43° 14304. 

The standard lest gives a probability of p(x^) < 
10~^ that the co-rotation straight line fits the data points 
when we use all 17 objects. Here we have considered the 
errors of the rotational period only, the errors of Porb are 
supposed to be small (usually they are <2%). The straight 
line cannot be rejected as a fit to the data at more than 
90% confidence level, p(x^)=0.127, when the two deviating 
objects (CD-43° 14304 and RS Oph) are excluded. The 
statistic tests how well the data are described by the model. 
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assuming that all the deviations are due to measurement 
errors and not to intrinsic scatter. 

When intrinsic scatter is allowed, as we expect to be the 
more realistic situation, one can apply the weighted least 
squares estimator (e.g. Akritas & Bershadi 1996) to find the 
slope {(3) and the intercept (a) of a linear fit in the form 
Pro* = a + (iVorb to the data points with measurement er- 
rors and non-negligible intrinsic scatter. When all objects 
are included, we get a = 204 ± 49 and /? = 0.14 ± 0.07, 
which is not consistent with Prot=Por6 {a = Q, (5 = 1.0). 
The situation changes if CD-43° 14304 is removed from the 
sample, then a = -110 ± 77 and 13 = 0.93 ± 0.16. When 
the deviating recurrent nova RS Oph is removed as well, we 
obtain a = -32 ± 78 and /3 = 1.01 ± 0.17. The last result is 
fully consistent with Vorb=Prot line (a = 0, /3 = 1). 

Although the synchronization line seems to fit well the 
data, we still cannot rule out the possibility that one or 
more outlying objects may determine a false correlation be- 
tween the variables. Therefore, a rank correlation test ap- 
pears appropriate. The Spearman rank correlation coeffi- 
cient between the two periods is 0.61 when all 17 objects 
were included, and about 0.82 for 16 and 15 objects only. 
This result implies a significant correlation between the vari- 
ables, since the p-values are less than 0.02 and 0.002 re- 
spectively. The similar Kendall-tau test gives similar results, 
with p-values less than 0.01 for all cases. 

The objects that deviate significantly from the 
Por6=Prot liue are RS Oph (a peculiar SS as noted above) 
and CD-43° 14304 (with possible error in the v sini mea- 
surement, see SectO. Taking into account the uncertainties 
of their Prot, one finds that the probabilities for a random 
deviation are less than 0.01 for them (see Table [5]), while for 
all the other 15 stars this idea cannot be rejected at the 0.01 
significance level. In other words the null hypothesis that all 
S-type SSs with well measured v sin i are synchronized (ex- 
cluding RS Oph and CD-43° 14304) cannot be rejected at 
the 99% confidence level. 

An additional test that may give some clues about the 
extent of synchronization of the periods is the K-S test of 
how much (Prot— Porb)/o" deviates from the normal distri- 
bution. Even when all 17 stars are included, the p-value of 
the K-S statistics is slightly less than 10%. When only the 
above 15 objects are included (see Fig|3}, the distribution 
becomes much narrower (mean =0.15, cr=1.23, K-S statis- 
tics is 0.38), with a standard deviation comparable to the 
measurement errors of P^ot, < ojProt >~ 0.53. This means 
that the hypothesis that the sample comes from a normal 
distribution cannot be rejected statistically. 

We see that the null hypothesis for synchronization of 
the red giants in symbiotic stars cannot be rejected statisti- 
cally (except probably for RS Oph). 




Figure 3. The distribution of (Prot-Porb)/'^ for 15 symbiotic 
stars and the fitted gaussian (cr is the error of Prot). The K-S 
statistic gives only a 38% probability that the parent distribution 
deviates from the normal one (not statistically significant). 



Table 5. The deviations of the individual objects from the line 
Prot =Porb- In the second column the deviations are given in 
units of (T, where rr is the individual error of Prot- The third 
column gives the probability for random deviation. To reject the 
null hypothesis (synchronization) at 99% confidence level, p{x^) 
has to be smaller than 0.01. 



Object 


deviation 











SY Mus 


2.2 


0.026 


AG Dra 


1.4 


0.175 


BX Mon 


1.2 


0.223 


TX CVn 


0.5 


0.605 


T CrB 


1.7 


0.082 


V443 Her 


0.8 


0.423 


CI Cyg 


0.7 


0.489 


AG Peg 


0.6 


0.569 


V1329 Cyg 


0.5 


0.651 


BP Cyg 


0.9 


0.361 


AR Pav 


1.2 


0.233 


V343 Ser 


2.2 


0.028 


BD-21°3873 


0.8 


0.399 


CD-43° 14304 


15.9 




FG Ser 


0.5 


0.610 


RS Oph 


6.7 


<io-« 


RW Hya 


0.4 


0.655 



9 CLUES TO THE ORBITAL PERIODS 

Up to now, out of 188 SSs, the orbital elements and binary 
periods are well known for ~40 objects only (and they are 
all S-type) . The derived orbital periods are in the range 200 
- 2000 days (Mikolajewska 2003). 

Because the orbital periods of the majority of SSs are 
unknown, an indirect method to obtain Porb is to mea- 
sure i; sini. If the mass donors in SSs are co-rotating 



(Prot =Porb), wc cau find clues for the orbital periods via the 
simple relation Porb Vrot ~ 27rRj, where Porb is the orbital 
period, Vrot and Rg are the rotational velocity and radius of 
the giant, respectively. The underlying suppositions are: (1) 
corotation (see Sect. [5]) and (2) rotational axis of the red 
giant perpendicular to the orbital plane (see Sect. IT2l and 
references therein). 

It could be useful in the case of eclipsing binaries, 
where sini ~ 1 and v sini~Vrot, however in such cases it 
is easy to find Porb with photometry. If the inclination is 
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unknown, we can only put an upper limit, Port <Puh where 
Pui ~ 27rRg/w sini. To see how useful this predictor is, we 
calculated the upper limits for the objects with known or- 
bital periods, and we get Pui/Porb ~ 0.5 — 1.5. 

For the objects with unknown periods, the upper limits 
for Port calculated in this way are given in the last column 
of Table (2] These upper limits are in the interval 100-1500 
days. Most of them are similar to those of the measured 
orbital periods in S-type SSs. However, it seems that for AS 
255 and Hen 3-1674, Porb could be as short as ;^150 days. 



10 CONCLUSIONS 

We have observed 30 S-type symbiotic stars with the FEROS 
spectrograph. We have measured the rotational velocities of 
the mass donors for 29 of them by the means of the CCF 
method (for V3804 Sgr we did not get a meaningful CCF). 
The results of the CCF method have been checked by the 
FWHM measurements. The main results are: 

(1) The projected, rotational velocities of the cool com- 
ponents in S-type symbiotic stars are from 3.5 kms~^ up 
to 52 kms~^ and 90% of them are in the interval from 4.5 
kms~^ to 11.7 kms~^. 

(2) In our sample of 17 S-type SSs with known orbital 
periods, 9 are synchronized within the measurement errors 
(1-a level). If we exclude the doubtful object CD-43° 14304 
and the recurrent nova RS Oph, all remaining 15 objects 
are synchronized within the 3-a level. In other words, 
the null hypothesis that these 15 objects are synchronized 
can not be rejected statistically at the 99% confidence level. 

(3) Among all objects with v sini measured by us the 
deviation from synchronization is statistically significant (at 
the 99% confidence level) only for RS Oph. The red giant in 
this peculiar object seems to rotate faster than the orbital 
period. 

(4) For 22 S-type SSs with unknown parameters, we 
give clues as to what their orbital periods could be. 

In future it will be interesting: (i) to measure the pro- 
jected rotational velocity of the cool giants in more symbiotic 
stars; (ii) to compare their rotational velocity with that of 
the isolated giants and those in other binary systems. 
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APPENDIX 

Here we give a few examples of our spectra (Fig. U) and 
graphic representation of the CCF for all SSs observed in 
this paper (Fig[5|. 
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Figure 4. A few examples of our spectra, together with the mask in the interval AA6400 — 6500 A(upper panel), and with a synthetic 
spectrum in the interval AA8740 — 8850 A. The spectra are plotted in increasing order of rotation: V4074 Sgr (M4III, v sini=3.5 kms~^), 
CD-43°14304 (K5III, v sini=7.2 kms"!), V919 Sgr (M2III, v sini=9.1 kms'^), RS Oph (MOIII, v sini=11.7 kms"!). The synthetic 
spectrum corresponds to MOIII and v sini=5 kms"'^. 
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Figure 5. CCF using KO numerical mask. In each panel are plotted the relative intensity of CCF (heavy line) and the fit versus radial 
velocity for the SSs observed in this paper. The measured widths of the CCF are given in Table |2] 
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Figure 5. continuation 



